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Abstract 

We have measured the neutron electric dipole moment using spin rotation in 
a non-centrosymmetric crystal. Our result is d n = (2.5 ± 6.5 stat ± 5.5 syst ) • 
10~ 24 ecm. The dominating contribution to the systematic uncertainty is sta- 
tistical in nature and will reduce with improved statistics. The statistical sensi- 
tivity can be increased to 2- 10 -26 ecm in 100 days data taking with an improved 
setup. We state technical requirements for a systematic uncertainty at the same 
level. 
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1. Introduction 

Electric dipole moments (EDMs) of elementary particles belong to the most 
sensitive probes for CP violation beyond the Standard Model of Particle Physics 
[1]. Constraining or detecting EDMs of different systems allows to gather ex- 
perimental information about models for new physics that is complementary to 
high energy physics data. 

For the neutron EDM (nEDM), the most sensitive results 0, Q were ob- 
tained using ultracold neutrons and Ramsey's resonance method. See Q for 
a recent review of measurements using free neutrons. Measurements using the 
interaction of neutrons with the atomic electric field in absorbing matter were pi- 
oneered by Shull and Nathans [H . Abov and colleagues @ first discussed a spin- 
dependent term in the scattering amplitude for neutrons in non-centrosymmetric 
non-absorptive crystals. This term is caused by the interference of nuclear and 
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spin-orbit structure amplitudes. Forte Q proposed to search for nEDM related 
spin rotation in non-centrosymmetric crystals due to such interference. In Q it 
was pointed out that such interference leads to a constant interplanar electric 
field affecting the neutron during all time of its passage through the crystal. 
This field was measured first in |9j. The interference of nuclear and spin-orbit 
structure amplitudes was tested by Forte and Zeyen [Io[ by spin rotation in a 
non-centrosymmetric crystal. 

Here we present the first measurement of the nEDM based on an improved 
version of this method. Preliminary results and a detailed description of our 
method with references to earlier work have been published in conference pro- 
ceedings (ll. 12 1 . 



The statistical sensitivity of any experiment to measure the nEDM is deter- 
mined by the product Et^/N , where E is the value of the electric field, r the 
duration of the neutron's interaction with the field and N the number of counted 
neutrons. New projects to measure the nEDM with UCNs aim to increase the 
UCN density and thus N by orders of magnitude and exploit in some cases the 
higher electric field obtainable in liquid helium compared to vacuum HHll. In 



contrast, experiments with non-centrosymmctric crystals exploit the interplanar 
electric field. For quartz, this field was measured to be E m 2 • 10 8 V/cm Q, 
several orders of magnitude higher than the electric field achievable in vacuum 
or liquid helium. Furthermore, the statistical sensitivity of the method profits 
from the higher flux of the used cold neutrons, compared to UCNs available 
today. These factors compensate the shorter interaction time of the neutrons 
with the electric field, ultimately limited by the absorption in the crystal. In 
[l5l | we have demonstrated that a statistical sensitivity of ~ 6 • 10~ 25 ecm/day 
can be obtained, comparable to the most sensitive published UCN experiments 
@,|. 

The results in 1 1 51 ] were obtained in Laue geometry that is limited by s\ 
tematics [l6j]. Here we exploit the Bragg geometry first proposed by Forte 



See [11| for a detailed comparison of the two schemes. 



2. Method 

We consider spin rotation for neutrons close to the Bragg condition for the 
crystallographic plane g in a non-centrosymmetric crystal. These neutrons are 
exposed to an electric field E = E g ■ a where E g is the interplanar electric field 
for plane g and a describes the deviation of the neutron from the exact Bragg 



condition (see for details). A nonzero nEDM d n results in neutron spin 



rotation by the angle 

2Ed n L 

</?EDM = — 7 , (1) 

nv± 

where L is the length of the crystal and v± the component of the neutron veloc- 
ity perpendicular to the crystallographic plane (L/v± is the time the neutron 
interacts with the field). By changing the deviation from the exact Bragg con- 
dition, a, value and even sign of the effective electric field and the resulting spin 
rotation </5edm can be selected. 
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On the other hand, the electric field causes a Schwinger magnetic field in the 
rest frame of the neutron, Hs = [E x v]/c, resulting in the spin rotation angle 



_ 2fi n H s L _ 2E f i n Lv\ l ^ 

hv± chv± ' 

where /Lt n is the neutron magnetic moment and the neutron velocity compo- 
nent parallel to the crystallographic plane. 

As E and H§ are perpendicular to each other, <^edm and ips can be separated 
by three-dimensional (3D) polarisation analysis. Furthermore, tps vanishes for 
Bragg angles of n/2 (v\\ = in Eq. ||2J)). This is used to suppress the effect 
due to Schwinger interaction: The crystal is aligned such that the interplanar 
electric field is parallel to the central neutron velocity, defining the Z direction 
of a coordinate system. The incident neutron polarisation is aligned in X (or V) 
direction. </5edm is measured in the X-Y plane. A residual Schwinger magnetic 
field (for neutron trajectories deviating from the Z direction or in case of a 
slight misalignment of the crystal) has its largest component in the X-Y plane, 
thus creating a polarisation component in Z direction. Thus, ip$ can be derived 
from the Z component of the outgoing polarisation vector. In summary, we 
measure <^edm from the component Pxy of the polarisation tensor and the 
residual Schwinger effect from the components Pxz and Pyz- Pyx, Pzx and 
Pzy serve for control purposes. 

In first order, the difference of the polarisation tensors for positive and neg- 
ative effective electric fields is: 



AP = g n r 



TO iJ^ivi, . To 



To 



EDM u - I n — -rn s 



(jjx At 
TO 



(3) 

where To = (r+ + t_)/2, At = (t + — t_)/2, and t + and t_ are the times 
the neutrons stay in the crystal for the positive and the negative electric field, 
respectively (the neutron velocity is slightly different for ±a). H l are the com- 
ponents of the residual magnetic field and H§ the components of the Schwinger 
magnetic field Hs- g n = 2^ n /h = 1.8 • 10 4 G _1 s _1 is the neutron gyromagnetic 
ratio, Aedm = Ed n / [i n the effective magnetic field corresponding to the electric 
field E. For E = 1 • 10 8 V/cm and d n = lCT 26 ecm, H EDM = 1.7 ■ 1CT 7 G. 



3. Experimental setup and procedure 

The experiment was carried out at the cold neutron beam facility PF1B 



17j of the Institut Laue-Langevin. A scheme of the setup is shown in Fig. [TJ 



Neutrons were wavelength-preselected by a pyrolytic graphite monochromator 
(adjusted to 4.91 A, resolution about 1%) and spin polarised to about 98% 
by a super mirror polariser. A resonance spin flipper permitted to flip the 
neutron polarisation. Cryopad (l8l . ITjj was used for 3D polarisation analysis. 
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Monitor Beamstop 




Figure 1: Scheme of the experiment. The neutrons come from the right. See Section l3l for 
details. 

The direction of the incident beam polarisation was oriented and the measured 
projection of the outgoing polarisation vector selected by nutators (in the X-Y 
plane) and precession coils (inside the Meissner screen, for the Z component). 
The polarisation of the transmitted beam was analysed by a 3 He spin filter cell 
(20j (measured polarisation value of the unperturbed beam between 75% and 
87%, depending on the 3 He cell). The currents in the precession coils were 
optimised experimentally for the wavelength of the used neutrons. 

We used the (110) reflection of a perfect quartz crystal of 14cm length. The 
effective angular mosaic spread of the crystal was u m ~ 1" . The neutrons with 
a well-defined deviation from the Bragg condition for this crystal were selected 
behind the 3 He spin filter cell by back-reflection at a second quartz crystal 
oriented parallel to the first one and at a slightly different temperature. A shift 
of the neutrons by one Bragg width (AAg/A w 10~ 5 for the (110) plane of 
quartz) corresponds to a temperature difference AT w ±1K (linear coefficient 
of thermal expansion for quartz £ = AL/L « 10 _5 /K). Note that the absolute 
temperature of the crystals is not important. The crystals were temperature 
stabilised in individual housings connected by a common water circuit. The 
temperature of the second crystal was varied by Peltier elements. 

The back-reflected neutrons were separated from the beam by a pyrolytic 
graphite monochromator (PG (002)) with a reflectivity of about 50% and di- 
rected to a position sensitive detector (nPSD). The nPSD allowed us to analyse 
data depending on the deviation of the neutron's path from the Z direction: 
Larger deviations correspond to larger uii and thus ipg. This monochromator 
unavoidably reflects about 50% of the neutrons after the spin filter, before the 
second quartz crystal. These neutrons were used to monitor intensity and po- 
larisation. 

During installation, the crystals were pre-aligned optically using an autocol- 
limator. Final alignment was performed with neutrons, allowing for the possibil- 
ity of slightly different lattice constants of the crystals at the same temperature: 
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The intensity of the neutron beam backscattered by the second crystal has a 
minimum if the Bragg condition of the second crystal (varied by the temperature 
difference) corresponds to the Bragg condition of the first one. If the crystal 
planes of both crystals are parallel, this minimum has the smallest width. The 
alignment with the minimum width was found by scanning the temperature dif- 
ference for different angles (in X and Y direction) of the second crystal. Data 
show that the lattice constants were equal for AT = — 0.8K [ll[, which can 
be attributed to slightly different impurities of the two crystals. From previous 
measurements [2lj we selected the working points with maximal electric field 
magnitude AT_ = -2.0 K and AT+ = +0.4 K. 

The different measured combinations of electric field direction and com- 
ponents of the polarisation tensor are Py±, where the index ± refers to the 
temperature difference of the crystals AT± . The indices i and j indicate the di- 
rections of the incident and outgoing neutron polarisation selected by Cryopad. 
Each component Pij± was measured for 40 s by flipping the incident neutron 
polarisation every 1-4 s with the resonance flipper (states j" off, J, on), in a 
drift-compensating scheme fiititti- The different states Pij± were measured 
in drift-compensating cycles, one for the EDM effect, Pxy-, Pxy+, Pxy+, 
Pxy-i Pyx-, Pyx+, Pyx+ ,Pyx-, and one for the Schwinger effect, Pxz-, 
Pyz-, Pxz+, Pyz+, Pyz+, Pxz+, Pyz-, Pxz-- EDM data were taken 2/3 
of the total time, Schwinger data due to the smaller required statistics only 1/3. 
The dominating drift originated from the depolarisation of the 3 He cells with 
time (time constants large compared to a single cycle). 

4. Results 

From the measured data we determined the elements of the difference of the 
polarisation matrices, AP^ = l/2(Py + — Pij-). The Schwinger data showed the 
expected linear variation over the nPSD: neglecting the spatial extension of the 
beam, the deviation of a neutron from the centre of the nPSD is proportional 
to vu. Fig. [2] shows a linear fit to the data APzy ■ This allows us to derive the 
electric field acting on the neutrons: 

P cxp = (0.7 ± 0.1) • 10 8 V/cm. (4) 

This value is consistent with our preliminary result obtained for the Bragg angle 

e B = 86° [23. 

Fig. [3] (top) shows an example for the spatial dependence of the matrix 
element APyx sensitive to the nEDM. The summary of all accumulated data 
for the nEDM is presented in Fig. [3] (bottom). The average value for the nEDM 
spin rotation is 

^ EDM = (0.9 ±2.3) ■ 10~ 4 rad, (5) 

corresponding to the nEDM (sec Eq. flTJ) d n = (2.5 ± 6.5) • 10~ 24 ecm. 

The slope of the spatial dependencies of APyx and APxy is related to 
the Schwinger effect. It was found to be consistent with zero and used to 
estimate a related systematic uncertainty of ad a> s < 5 • 10~ 24 ecm. Note that 
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Figure 2: Measured dependence of APzy on the spatial coordinate Y of the nPSD. 

this uncertainty is smaller than the statistical uncertainty of the nEDM taken 
from the same data and would reduce with better statistics. It includes the 
error caused by imperfect angular accuracy of the 3D polarisation analysis (if the 
analysed Z component is not perpendicular to the X-Y plane, the corresponding 
projection of the Schwinger effect would contribute to APxy and APyx)- 

The residual magnetic field inside Cryopad was H ~ (1—2) mG. The reversal 
of the electric field leads to a small change of the neutron wavelength and thus to 
a difference of the time the neutron stays inside Cryopad. This time difference 
was measured directly: At/t < 5 • 1CP 4 . The combination of the residual 
magnetic field and the time difference results in a systematic uncertainty of 
ad nt H < 4.5 ■ 1CT 25 ecm. 

The time stability of the experiment was controlled by the beam monitor, 
which registered neutrons behind Cryopad without interaction with the second 
quartz crystal. The average electric field acting on these neutrons was about 
zero (taking into account the wavelength distribution in the neutron beam, it 
was ~ 10 3 times less than the field for the neutrons that would be selected 
by the second quartz crystal). The difference in polarisation measured by the 
beam monitor between the two working points, AP < 2 ■ 1CP 5 , can be used to 
limit a possible systematic effect due to beam fluctuations in time, (jd n ^ lm c < 
1(T 24 ecm. 

We have observed a spatial shift of the neutron beam at the nPSD corre- 
lated with the working point AT± (see Fig. 2k), caused by an imperfect align- 
ment of the two quartz crystals, and a spatial dependence of the polarisation 
(see Fig. HJd), probably due to the cylindrical shape of the outgoing window 
of Cryopad and the presence of a residual magnetic field. The combination of 
both effects causes a systematic error of AP S = 1/2 AY S SP/6Y S . For the fi- 
nal crystal alignment AY S = 0.03 cm this corresponds to AP S = 4 ■ 10 -5 and 
ad n , s = 2 ■ 10~ 24 ecm. 

Spin rotation caused by the weak interaction cancels for the two crystal 
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Figure 3: Top: Measured dependences of APyx on the spatial coordinates X and Y of the 
nPSD. Bottom: Results of the different measurements of APxy and APyx, used to derive 
the nEDM. 



temperatures in first order. The cancellation is incomplete because of a cross 
term of weak and Schwinger rotation, but the uncompensated part is below 
10 -4 of the weak spin rotation and can be neglected. 

Summing these independent systematic errors we obtain the final result 

d n = (2.5 ± 6.5 stat ± 5.5 syst ) ■ 10" 24 ecm. (6) 

Note that the systematic error is mainly statistical in nature. This result is 
about two orders of magnitude less precise than present experiments using UCNs 
and Ramsey's resonance method. 



5. Outlook 



The statistical sensitivity of the experiment described above was 1.6 • 
10~ 23 ecm/day. In (ll| we have shown that the sensitivity can be improved by 
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Figure 4: Spatial distribution of the registered neutron intensity a) and measured polarisation 
b) for the two working points AT± . 



Table 1: Experimental parameters for a systematic uncertainty of ad n ~ 10 26 ecm. 



Parameter 


Value 


Interplanar distance variation Ad/d 


C lQ- b 


Residual magnetic field 






Value 


- 1(T 4 G 


Time stability 


- 10-°G/h 


Accuracy of 3D polarisation analysis 


- Kr a rad 


Angular crystals alignment 


- 0.01° 


Precision of temperature control 


- 0.01K 


Flatness of Cryopad windows 


- l(T 4 rad 



a factor 65 to about 2.5- 10~ 25 ecm/day, essentially by increasing the divergence 
acceptance of the installation, the beam size, and the crystal length. This im- 
proved sensitivity is comparable to state-of-the-art nEDM measurements using 
UCNs and permits a statistical uncertainty at the 10 -26 ecm level in a beam 
time of 100 days. A further improvement may be possible by using crystals 
with higher interplanar electric fields. However, these crystals are presently not 
available in the required quality and size. 

In Table [1] we list the experimental parameters that are required for an ade- 
quate systematic uncertainty. The values are achievable. Note that a variation 
of the interplanar distance in the crystal of the order of the Bragg width or 
larger does not simulate an nEDM but reduces the magnitude of the effective 
electric field. 

We conclude that neutron spin rotation in quartz in Bragg geometry per- 
mits to measure the nEDM with a precision comparable to existing UCN ex- 
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periments but using a completely different method. Efforts are now being made 
to produce larger crystals and a version of Cryopad optimised for that kind of 
measurements. 
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